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ABSTRACT 

Comparing the frequency of typical events with that of unusual events allows 
one to test whether the cosmological density distribution function is consistent 
with the normally made assumption of gaussianity. To this end, we compare the 
consistency of the tail-inferred (from clusters) and measured values (from large- 
scale flows) of the rms level of mass fluctuations for two distribution functions: 
a Gaussian, and a texture (positively-skewed) PDF. We find that if we average 
the recent large-scale flow measurements, observations of the rms and the tail at 
the 10/i _1 Mpc scale disfavor a texture PDF at ~ 1.5a in all cases. If we take 
the most recent measurement of the rms, that from Willick et al. (1997b), the 
comparison disfavors textures for low Qq = 0.3, and disfavors Gaussian models if 
Qo = 1 (again at ~ 1.5a). Predictions for evolution of high temperature clusters 
can also be made for the models considered, and, as is known (e.g., Henry 1997), 
strongly disfavor Qq = 1 Gaussian models, while we find Qq — 1 marginally 
disfavored in texture models. Taking the suite of tests as a whole, and using all 
of the quoted data, it appears that textures are strongly disfavored and only the 
low Qq Gaussian models are consistent with all the data considered. But given 
evidence for the internal inconsistency of the observational data, had we only 
used the recent Willick et al. results, the strength of our conclusion would have 
been reduced to the ~ la level. 

Subject headings: cosmology: large-scale structure of universe — galaxies: clus- 
ters and statistics — methods: analytical 
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1. Introduction 

All current scenarios for the formation of structure in the Universe assume that cosmic 
structures grow from initially small seed perturbations. A fundamental assumption underly- 
ing a large class of these theories is that the initial Fourier modes have random phases. This 
implies that the statistics of the initial density field pj(x) are fully specified by the correlation 
function £$(r) = (pj(x)pj(x + r)) x , or equivalently its Fourier transform, the power spectrum 
Pi(k). For any given smoothing length, then, the one-point probability distribution function 
(PDF) of the initial density field pi is Gaussian. A measurement of the frequency of rare 
events (clusters of galaxies) coupled with a measurement of the level of rms fluctuations can 
tell us directly whether or not the distribution of density fluctuations satisfy the Gaussian 
hypothesis. 

While theories based on the inflationary paradigm generally predict Gaussian initial 
conditions, there are many alternatives which do not. Global textures and other topological 
defects, for example, predict a strongly non-Gaussian distribution for the initial perturba- 
tions (e.g., |Gooding, Spergel, fc Turok 1991| ; Park, Spergel, fc Turok 1991| ; |Gooding et al.| 



1992| ; Pen, Spergel, fc Turok 1994] ; Pen et al. 1997] ). Global textures arise when a global 



non-Abelian symmetry, such as SU(2), is spontaneously and completely broken. When these 
defects collapse, energy gradients in the Higgs field perturb the metric and induce gravi- 
tational fluctuations in the matter fields. These texture knots, then, act as "causal seeds" 
of structure formation. We take the texture model as our standard non-Gaussian model, 
since it is well studied and still viable as a theory of structure formation. In general, texture 
and Gaussian models give similar predictions for the shape of the power spectrum and the 
level of CMB fluctuations ( Pen, Spergel, fc Turok T994] ; Pen fc Spergel 1995Q . In fact, it is 
often noted (e.g., Park, Spergel, fc Turok 19911) that texture and Gaussian CDM models are 



most different in their PDFs. If it could be demonstrated that the initial density field were 
Gaussian on some suitable scale, the texture model would be directly falsified. 

The standard method for testing Gaussianity relies on measuring the iV-th order reduced 
moments, or cumulants A at of the present-day PDF, such as skewness or kurtosis, from 
redshift or peculiar velocity surveys (cf., Strauss & Willick 1995 for a review). The variance 
is a 2 = A 2 = (S 2 ), the skewness is A 3 = (<5 3 ), and the kurtosis is A 4 = (<5 4 ) — 3cr 4 , where S 
is the density contrast p/(p) — 1. For a Gaussian random field, all the cumulants of order 
N > 2 are equal to zero. The growth of Gaussian initial perturbations under gravitational 
instability naturally produces non-Gaussianity in the distribution function of the density 
field. In particular, the cumulants follow a scaling relation Xn = StvA^ -1 , where the Sn 
are independent of scale in the mildly non-linear regime for power-law power spectra (Fry 
1984ab, Bernardeau 1992), and can be calculated from (N — l)-th order perturbation theory 
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(Bernardeau 1994). Significant departures from scale-invariance are taken as signs of non- 
Gaussianity. Unfortunately, the higher-order moments are very sensitive to finite-volume 
effects and non-Gaussian noise in the data (c.f. the discussion in [Kim fc Strauss 1997] and 
references therein). 

Since many non-Gaussian distribution functions have long tails, it is better to directly 
compare the rms and the tail of the PDF, rather than relying on the moments, which are 
integrals over the PDF. Given a distribution function, the level of tail fluctuations maps 
directly to a level of rms fluctuations; the tail-inferred rms, then, must be consistent with 
the directly measured rms. Furthermore, one should look for measures which are, as closely 
as possible, indicators of the initial, linear density field. 

In this paper, we use the observed temperature distribution of clusters of galaxies as 
a measure of the tail of the distribution of initial density fluctuations, and galaxy peculiar 
velocities and redshift-space distortions as measures of the rms. Clusters of galaxies are the 
most massive virialized objects in the universe, with densities a hundred or more times that 
of the mean density of the universe, containing only a few percent of all galaxies. They are 
thus extreme excursions from the mean, and hence are measures of the high potential tail 
of the distribution of fluctuations. We compare this with the rms level of peculiar velocity 
fluctuations. The rms is not as sensitive to non-linear growth as are higher order statis- 
tics (Dekel, private communication). Both cluster and peculiar velocity observations are 
at roughly the same spatial scale, ~ 10/i -1 Mpc. Furthermore, they both measure gravita- 
tional potential fluctuations, and therefore have a similar dependence on Qq, which can be 
approximated by defining 

VR = ^O 6 , (1) 

where <7r is the rms mass fluctuation on tophat filter scale R. The remaining Q dependence 
is < 5% for the peculiar velocity and redshift distortion measurements, and < 25% in the 
range 0.3 < Qq < 1.0 for the cluster measurements. This is the quantity we will attempt to 
determine throughout our analysis, comparing estimates of it found by various methods. 

In §2, we examine various recent derivations of 7710 from the Mark III peculiar velocity 
catalog (Willick et al. 1995, 1996, 1997a), and from IRAS galaxy redshift space distortions. 
In §3, we infer rjio from the redshift-zero X-ray cluster abundances, assuming a Gaussian and 
texture PDF, using a modified Press-Schechter approach. We calculate the rjio separately for 
Qq = 1.0, cosmologically flat (A-dominated) Qq = 0.3, and open Q = 0.3. In §4, we discuss 
the consistency of the rms and the tail-inferred measures of 77 by performing a likelihood 
analysis of the Gaussian and texture hypotheses. We also discuss further constraints from 
observations of clusters at moderate redshift. We summarize in 55. 
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2. RMS Mass Fluctuations from Peculiar Velocities and Redshift Distortions 

It has been noted for more than a decade that the rms optical galaxy counts in 8/i _1 Mpc 
spheres is ((5N/N) 2 S ) ^ = a °P tical « l (p a vis fc Peebles 1983|) . But given that galaxy for- 



mation is still not very well understood, it has become customary to relate this measured 
variance to the more fundamental rms fluctuations in the matter density Ug by a bias param- 
eter b. For optical galaxies, fr°P tlcal is defined by <j° ptlcal = &°P tlcal or 8 . The bias, however, is 
clearly dependent on the galaxy type (for example, as we will see below al RAS = 0.69 ±0.04, 
[Fisher et al. 1994a|) , and may be dependent on scale. We wish to estimate the rms value of 
fluctuations independent of bias. To do this, we look at recent analyses of peculiar velocity 
and redshift survey data. 



2.1. Recent Measurements of RMS Fluctuations 

The peculiar velocity v is the deviation of the proper velocity of a galaxy from the local 
Hubble flow. Assuming that this measured velocity follows from gravitational potential 
fluctuations, we obtain from linear perturbation theory ( Peebles 1993| ): 

V • v = -HJ(n )S « -H o n° 6 5. (2) 

Here 5 is the mass density contrast S(x) = pii n (x)/p — 1, and the approximation is good to 
5%. Now consider filtering the velocity and density fields with a tophat of radius R, yielding 
the filtered fields vr and Sr. Then V • vr, the divergence of the filtered velocity field, is 
proportional to the filtered mass density field Sr. Taking the rms values of V • vr and Sr 
gives: 

^((v ■ v R ) 2 > Va » n^(sl) Va = ee U r. (3) 

The peculiar velocity field thus provides a direct measurement of mass fluctuations, appearing 
with the cosmological density parameter Q . Note that in linear theory, where individual 
Fourier components evolve independently, it makes no difference in which order one filters 
and takes the divergence of the velocity field. 

Seljak & Bertschinger (1993) made an early attempt to measure rjR directly using the 
POTENT method (Bertschinger et al. 1990) to reconstruct the mass density and peculiar 
velocity fields from the Mark II catalog (Burstein 1989). They derived a value of = 
1.3 ±0.3. More recently, several workers have used the Mark III catalog of peculiar velocities 
(Willick et al. 1995, 1996, 1997a) to derive values of rj. 

Zaroubi et al. (1997) perform a likelihood analysis of the peculiar velocities given families 
of CDM-type models. Because we are making a comparison only with clusters of galaxies, 
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which are measures of fluctuations on scales of roughly 10 h~ 1 Mpc, we only consider their 
COBE-independent results. Zaroubi et al. (1997) work with the measured peculiar velocities 
of individual galaxies, without any smoothing, so it is unclear to which scale they are most 
sensitive. We simply integrate their power spectra with 10/i _1 Mpc tophat smoothing to 
obtain the rms fluctuation on that scale. Their "maximum likelihood" model has 7710 = 0.70, 
and the models on their "65%-confidence" level contour span the range r] w = 0.61 — 0.78. 
We therefore take their measurement to be 7710 = 0.70 ± 0.09. 

Kolatt & Dekel (1997) use the POTENT technique (Dekel et al. 1990) to determine the 
quantity V • v# from the Mark III peculiar velocities, to compute the mass power spectrum 
directly, subtracting out a model for the noise power spectrum derived from mock cata- 
logs. Since POTENT applies a Gaussian filter of filter length 12/z _1 Mpc, equivalent to a 
~ 23/i -1 Mpc tophat filter (see Appendix), we integrate the functional form they give for the 
observed power spectrum, taking into account the uncertainties in the parameterization, to 
derive a value of rj 2 3 = 0.32 ±0.11. 

If we assume linear biasing, 5 ga iaxies = &<Wk matter, comparison of peculiar velocity and 
galaxy density field data via equation (|2]) allows one to determine the quantity ft = ^q 6 /6, 
the proportionality constant between the peculiar velocity and galaxy density fields for a 
particular set of galaxies. Note that if bias is independent of scale, ft should be as well. 
Following equation (|3|) and the definition of b, we find that for a comparison with IRAS 
galaxies: 

Vr — Piras&r,iras- (4) 
For ctrjras, we use the power law analytic approximation of Fisher et al. (1994a): 

720 IBV 1 

; r = SJel^rtlpc; 7 = 1-661^. (5) 



or 



27(3- 7 ) (4- T ) (6 -7) 



Several groups have determined ft from comparisons of IRAS density and peculiar ve- 
locity fields. Willick et al. (1997b) perform a likelihood analysis of Tully-Fisher observables 
of the Mark III data given a prior velocity model, and find ft mAS = 0.49 ± 0.07. Their 
analysis uses an effective Gaussian smoothing of ~ 4/i~ 1 Mpc and so corresponds roughly to 
an 8/i~ 1 Mpc tophat. Using cTg RAS = 0.69 ± 0.04, we obtain rj 8 = 0.34 ± 0.05. In another 
approach, Sigad et al. (1997) compare the POTENT-derived peculiar velocity divergence 
field to the IRAS galaxy density field (with effective tophat smoothing of ~ 23/i _1 Mpc due 
to the POTENT method) and derive ft = 0.89 ± 0.12. Multiplying by c^f AS = 0.29 ± 0.02 
leads to r/ 23 = 0.26 ± 0.04. 

Finally, as first pointed out by Kaiser (1987), one can use the anisotropy of the redshift- 
space correlation function of a redshift survey to infer the value of ft. Fisher et al. (1994b) 
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measured the anisotropy of the IRAS 1.2Jy sample at an effective scale of 10 — 15/t _1 Mpc, 
obtaining (3 = 0.45lo^g. Using the same procedure as above to convert (3 to rj, we obtain 
7712.5 = 0.21 ± 0.11. Cole et al. (1994, 1995) use instead the power spectrum anisotropy, 
at an effective scale of ~ 35/i _1 Mpc, to obtain (3 = 0.52 ± 0.13, resulting in r/35 = 0.11 ± 
0.03. Yet another method involves decomposing the redshift-space distribution into spherical 
harmonics. Fisher, Scharf, & Lahav (1994) (c.f. [Heavens fc Taylor 1995| ) show that the 



inferred value of (3 depends sensitively on the assumed power spectrum. Scaling their results 
for each of the power spectra they use, we find that the results for 77 roughly converge at 
SO/i^Mpc, with a value r/ 30 = 0.20 ± 0.03. 

Table [I] summarizes the values of r/# cited above, all using rms fluctuations, while Fig- 
ure ^ plots all the recent measurements, with their respective errors, at their effective tophat 
filter scales. Most current theories of structure formation predict that or is approximately 
a power law in the range 8 < R < 35/i -1 Mpc. We thus approximate T)r as a power law 
T)r oc R~ a , to allow us to extrapolate the results of Table [l] and Figure [l] to a common scale. 
In linear theory, a standard CDM model with n = 1 and Qq = 1 has a logarithmic slope at 
lO/i^Mpc of a « 1.0, while a tilted n = 0.7, low density fi = 0.3 CDM model has a « 0.6. 
An fl = 1 texture model (using the power spectra of Pen, Spergel, fc Turok 1994] and Pen| 
fc Spergel 1995| ) gives a ~ 1.2 at 10/i _1 Mpc. If galaxy bias does not change appreciably 



in this range, then a is the same as that for galaxies. Equation (|5]) implies an exponent 
a = 0.83 ± 0.12 for galaxy fluctuations. If the true value of a for mass fluctuations is much 
less (greater) than 0.83, then bias must decrease (increase) with increasing scale. When ex- 
trapolating the above results to different scales, then, we take three different power laws to 
reflect this uncertainty in the slope of the fluctuation spectrum, t]r cx R~ a : a = 0.83, 0.60, 
and 1.10. Table H] contains extrapolations of to 10/i _1 Mpc for each of these assumed values 
for a. It is striking that the different methods give answers which differ from one another by 
significantly more than the quoted errors — a sure sign that systematic errors of unknown 
origin are important. We now present a method for combining these measurements. 



2.2. An Error Model for the Different Measurements of i] 

Let us assume that each of the seven measurements (labeled i — 1 . . . 7) listed in Table [I] 
has an unknown systematic error /ij from the true value of rj. That is, each one actually 
measures rj + fii, rather than r/, where the stated statistical error q is the error in the 
measurement of 77 + /ij. We assume that the systematic error is drawn from a Gaussian 
distribution with zero mean and unknown variance 8 2 , and that the statistical errors q are 
also Gaussian distributed. The probability of a measurement x i: given unknown 77 and 9, is 
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given by 



V(xi\r],6) = J dniV{xi\Hi : ri)V{ni\6) 

(Xi - 7?) 2 



exp 



+ 3. 



(6) 
(7) 



2tt(0 2 + e?) 

The systematic error simply adds an error # in quadrature to each measurement error e^. 
If we define the likelihood function of the unknowns 77 and 9 in the usual way, we obtain 



In C(rj, 
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2tt( 



(9) 



^ 2 + 6? 

The maximum likelihood value for 77 is obtained by setting the derivative of equation (H) to 
zero: 

Xi \ 1 



-1 



V 



E 



# 2 + 



X 



E 



# 2 + e 2 



(10) 



Note that equation (|I0|) gives the traditional weighted average for 9 -C e», and an unweighted 
average for # 3> 6j. Since we are primarily concerned with 77, we numerically integrate the 
likelihood £(77, 9) over # to obtain the marginal likelihood over 77: 



^rms yl) 



C(r),9)d9. 



(11) 



In Figure [| we show the likelihood C vms (r]), each normalized to j drj C = 1, in the case of 
a =0.60, 0.83, and 1.10, along with a Gaussian distribution of the same 68% confidence level 
in the case of a = 0.83. Note the significant tail to the distribution relative to a Gaussian, 
reflecting the fact that the measurements are not consistent their individual normal errors. 
The values of rj W corresponding to the maximum of £ rms (i;), equation fllTD , along with la 
(68% confidence) errors, are given in the last row of Table |], for various a. For a = 0.83, 
the rms measurements of the velocity field give 77 10 = 0.44 ± 0.08 (ler). 

This model is only a rough attempt at quantifying the degree of systematic errors present 
in these measurement. The assumptions that each systematic error ^ is independent and 
normally distributed are certainly not completely correct. All of the analyses use the IRAS 
galaxy redshift catalog, and many of them use the same Mark III catalog, and hence their 
results must be somewhat correlated. Furthermore, systematic errors are rarely normally 
distributed. Perhaps one of the measurements is correct, and all the others are erroneous. In 
comparing rms and tail-inferred values of 77, we therefore also use the Willick et al. (1997b) 
result alone, which is the most recent result which explicitly filters at a scale close to the 
~ 10/i -1 Mpc scale from which clusters form. 
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3. Tail Fluctuations: Abundances of X-ray Clusters 

Galaxy statistics and peculiar velocities measure the rms fluctuations of the density 
field — i.e., fluctuations within ~ la from the mean. Most probability distributions look 
similar to Gaussians at this fluctuation level — they are peaked and convex near the mean. 
The more dramatic differences between Gaussian and non-Gaussian distributions come in 
their tails. Rich clusters of galaxies, the most massive virialized objects in the universe, 
offer the best measure of the tail of the probability distribution of the density field on scales 
of ~ 10/i _1 Mpc. Approximately 5% of the galaxies are within one Abell radius of a 
rich cluster ( [Bahcall 1996| ), and about 10% of the baryons reside in clusters (taking a BBN 
value of Qb = 0.0125/i -2 from Walker et al. 1991| ). Thus in the Gaussian scenario, they 



constitute ^ 1.5<r perturbations. They are well-observed by a number of techniques; in 
particular, their abundance as a function of X-ray temperature is well-measured (Henry & 
Arnaud 1991). Thus they provide an accurate estimate of the integral over the tail of the 
distribution function at the ^ 1.5cr level. 

The temperature function of rich clusters can be predicted from the statistics of peaks in 
a density field; it depends on the rms level of fluctuations or on the scale on which clusters 
form, as well as Qq, through the relation between comoving radius and mass. Comparing 
these predictions with observations constrains the quantity <Jr£Iq, v = 0.4 ~ 0.6 (cf., White, 
Efstathiou, & Frenk 1993; Eke, Cole, & Frenk 1996; Pen 1997, and references therein). The 
mass of a rich cluster within the 1.5/i~ 1 Mpc Abell radius is 10 14 — 2 x 10 15 /i~ 1 Mq, correspond- 

ing to a region with initial comoving radius R = [3M/(47rp G )] h = (4.4-12)fi ^/i^Mpc. 
Numerical simulations have confirmed that the cosmic abundance of these objects is deter- 
mined by the power spectrum at k~ x = 5 — 20/i _1 Mpc, so that their frequency provides a 
direct measure of or in this range (e.g., Bahcall & Cen 1992; |Fvrard et al. 1996] ; Cen 1997, 



in preparation). Hence the cluster abundances constrain r/io with only a weak dependence 
on Q . In §3.1 we derive the relation between the distribution of clusters as a function of 
temperature and the assumed density PDF. This requires a relation between the observed 
temperature of a cluster and the initial comoving radius from which it formed; this is given 
in §3.2. This relation depends on the formation epoch of the cluster, as derived in §3.3. We 
compare these results with the observed temperature distribution of clusters in §3.4 to derive 
the quantity r] 10 . 
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3.1. Press-Schechter and the Temperature Function 

We use the Press-Schechter (1974) ansatz (hereafter PS) to derive the number density 
of clusters. Many have done this in the Gaussian case (e.g., White et al. 1993), and here 
we make a straight-forward generalization to non-Gaussian distributions. In PS, an object 
of mass > M forms when the linearly-extrapolated density contrast filtered on a mass scale 
M, corresponding to a tophat comoving length scale R, exceeds the threshold 5 C ~ 1.69. 
The exact value of 5 C depends on flo and A, and can be derived assuming uniform spherical 
collapse (see [Eke, Cole, fc Frenk 1996|) , but varies by only a few percent for Qq > 0.2. The 



differential abundance from PS, given here as the comoving number density at the present 
of objects with initial comoving radius in the interval (R, R + dR), can be written as 



riRdR 



AnR 3 



d f°° 
dR Js c /a R 



dR, (12) 



where we assume that the mass distribution function P(y) only depends on y = 5/cr R . This is 
true in Gaussian theories; in the texture model, the seeding of density perturbations is scale- 
invariant, and thus the distribution function must also be scale-invariant. The prefactor 
/ = [/ °° P(y)dy}^ 1 normalizes P so that all the mass in the universe is accounted for. 
Assuming, as before, a power law dependence or oc R~ a gives 

3 a V f r>( \ _ S c , 1q n 
nR= 4^1l fP{y) > y = ~ R (13) 

The initial comoving radius is of course not an observable, so we will convert to the abundance 
as a function of present-day temperature: 

dT 3 ay , . t . 

n T — = --fP(y). (14) 

dR AttR 3 R j Kyj K J 



The observed X-ray temperature function gives the nr of equation ( |T4] ) as a function of 
temperature. In order to find y (and hence a and rj), we need the appropriate PDF for each 
model, and the relationship between the initial comoving radius and the temperature of the 
clusters (next section). 

In Gaussian theories, the distribution function is 

/P(y) = -LexpH/ 2 /2), (15) 

recalling that the factor / = [J °° P^dy]^ 1 = 2. For textures, we take the distribution 
function from the numerical simulations of Park, Spergel, & Turok (1991), which are the 



-10- 



only published simulations which are at scales appropriate for clusters of galaxies. We find 
that the high-?/ tail of their distribution gives: 

/P(y) = 2.2exp(-1.45y), y > 1.5. (16) 

Equation (|T6| ) is not valid for y < 1.5, but as we are only considering the tail of the distri- 
bution, it suffices for our purposes. We show the tails of the functions / P(y) for a Gaussian 
and texture PDF in Figure |3|. Note that the texture PDF has a significantly wider tail for 
y>2. 



3.2. The Radius- Temperature Relation 

We next consider the relationship between the initial comoving radius of the perturba- 
tion and the present temperature of a cluster. We make the standard assumption that a 
cluster forms by undergoing a spherical collapse to a singular isothermal sphere with a virial 
radius equal to half its maximum expansion radius (e.g., White et al. 1993, Eke et al. 1996, 
Pen 1997). The X-ray temperature simply reflects the depth of the gravitational potential 
well of the cluster. We assume a hydrogen mass fraction of 0.76. For clusters collapsing at 
redshift Zf, we obtain 

rp = ym P GM 



2r 



vir 




(19) 



lOh-m-pcJ \18tt 2 J \A C (0) 
= T (R)xcl>(z f ), (20) 

where the quantity in the square brackets is designated 0, and only depends on the back- 
ground cosmology and the formation redshift Zf. The function T Q (R) is the Radius- Temperature 
relation for collapse at Zf = 0, for a given background cosmology, thus <p(zf = 0) = 1. Equa- 
tion (|0|) was derived previously by Eke et al. (1996). The factor A c = 3M/(47rry ir p c ) is 
the ratio of the virialized physical mass density to the critical (not the mean) cosmological 
density at the time of collapse, and r vir is taken to be r max /2. In an Q = 1, matter dominated 
universe, A c = 187r 2 pa 178; for Q < 1, A c can be derived from uniform spherical collapse 
(see Eke et al. 1996). The quantity A c depends only on the background cosmology. Because 
A c and Q asymptote to 187r 2 and 1 respectively as Zf — > oo, and hence change by a factor 
< Qq 1 , the z/-dependence of <p is dominated by the factor 1 + Zf. 
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Evrard et al. (1996) have compared numerical simulations for Gaussian Q = 1.0 and 
fi — 0.2 models, and found that T scales as M 2//3 , as equation flUf ) implies. Fitting a 
power law to their results, and converting mass to comoving scale implies kT = (8.8 ± 
0.8 keV)(A c /187r 2 )°- 164 fio /3 (^/8/i _1 Mpc) 2 . This equation is consistent with equation (|TJ) 
for all three background cosmologies we consider, using the results of the next section to 
determine Zf and <f)(zf). The detailed simulations necessary to check the cluster temperatures 
for the texture scenario have not been done. Bartlett (1997) notes, however, that for any 
given density profile of a spherical collapse, the power law T oc M 2/3 oc R 2 should be the 
same — the only difference might be in the constant of proportionality (e.g. through a 
different ^-dependence) . 



Taking the derivative of equation (^0|) with respect to R gives 



dT rrf D \d(f> r ,dT T ( d\n(f)\ 

dR =To{R) dR + ^iR =R\ 2+ d^R)' (21) 

where the redshift of formation Zf must be calculated as a function of R. Finally, inverting 
equation (|2"0"|) gives R as a function of the temperature: 



R = 10/,-Mpc fio - 1 Af^ 7 ) 1/2 fi^V /8 x 



8.6 keV/ V A c(°) 



Ac( Z f) ^0 J (l + Zf ) l /l 



= i? o (T)x0[^(y)]- 1 /2 ) (23) 

where <fi is defined as above, the factor in front of the square brackets has been replaced by 
Rq(T), the Radius- Temperature relation at Zf = 0. Note that equation ( P^D so defined is 
only valid when Zf is known. 



3.3. Formation Epoch of Clusters in PS 

We now have another parameter to consider, the redshift of formation Zf. The standard 
assumption is to assume that the cluster has "just formed" at the redshift z at which we 
observe it, i.e., that Zf = zq. In the Qq = 1 Gaussian CDM scenario, this assumption is valid, 
since clusters form late and accretion continues to the present (e.g., Lacey & Cole 1993). 
If O < 1, however, linear growth of perturbations freezes out when the universe enters its 
phase of free expansion at 1+z ~ 1/Qq (see Peebles 1993 ). Therefore, a given structure must 



collapse at a higher redshift than in the Qq = 1 case in order to form. If there is a sufficiently 
positively-skewed non-Gaussian tail to the distribution function, as there is for textures, a 
larger fraction of the universe collapses non-linearly early on (see Bartlett et al. 1993 and 
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Pen et al. 1994). Furthermore, for textures and other "causal-seed" models, whose density 
perturbations are compensated, there is a limit to the final mass of a collapsing structure, 
defined by the size of the causal horizon at the time the perturbation is seeded. Because the 
integrated overdensity is zero on larger scales, no further accretion can occur (Bartlett et al. 
1993). 



We now derive the mean redshift of formation (z/), and use that redshift in equation (p0[) 
to determine the typical temperature of clusters as a function of R. The net differential num- 
ber density of clusters which form in redshift interval (z, z + dz) follows from equation ([13]) : 

where 

v(z) - W = S c (0) S c (z) 

= yo JmWY (26) 

Here dq refers to the value of y(z) at z = 0, D(z) is the linear growth factor normalized 
to unity at redshift zero, and S c (z) is the critical linear density of collapse at redshift z. 
For clusters observed at redshift zq, the mean redshift of formation (zf) ZQ is obtained by 
averaging the redshift z weighted by equation (|24T): 



\ Z fho = roo dnlj^ ( 27 ) 

h dz az 

- ^ + ^ikm £ yiz)p[v{z)]dz (28) 

where the second line comes from integrating the numerator by parts, integrating the de- 
nominator, and cancelling terms. The mean redshift of formation of clusters at z is a 
function only of the background cosmology through S c (z) and D(z), and the current value 
of yo through equation (p6|). We take the formation redshift required in equation (p0|) to be 
(zf). When using equation ( p8[) with the present epoch abundances below, we set z = 0. 



3.4. Determining 77 from X-ray Cluster Surveys 

Using equations (|20|) and (pl~D in equation (p!4|) gives 

m ( a dhi(t)\ 3 „ . . . . 



13 



where y is now understood to be its present day value, and we have used = 



a. 



because of our power-law assumption. Replacing R with equation ( p3]) and separating those 
quantities depending on T from those depending on y gives 



(30) 

1+ 2 d\ny 

We need now only to relate observations of clusters to the left-hand-side of equation (|30|) . 

There are two virtually complete X-ray flux-limited surveys of rich galaxy clusters at 
low redshift. The first was compiled by Henry & Arnaud (1991), and the second is the 
"X-ray Brightest Abell Cluster Survey" (XBACS) from the ROSAT all-sky survey ([Ebeling 



et al. 199"6|). The Henry & Arnaud (1991) survey is essentially contained within the XBACS 



survey. An unbiased estimator of the number density of clusters in the temperature interval 
(T, T + AT) for a flux-limited survey is given by 

n T AT= £ V~^, (31) 

T<Ti<T+AT 

where V mauX ^ is the maximum volume in which the cluster of temperature Tj (in keV) could 
be detected given the flux limit and geometric boundaries of the survey. In this case, we 
require a geometrical boundary at galactic latitude |fe| > 20°, and redshift z < 0.1. The 
full XBACS sample takes the maximum redshift to be z max = 0.2. The Poisson variance, 
following Pen (1997), is given by: 

a 2 (n T AT)= £ V^. (32) 

T<Ti<T+AT 

The cumulative temperature functions n(> T) for the two samples are shown in Figure |[ 
The two temperature functions are quite consistent in the overlapping temperature range 
> 3keV, so we use this part of the (10 x) larger XBACS survey for our analysis. 

Substituting equation ([31]) gives, finally, 

87tRl(T) 4>layfP(y) 




l " 1,1 11111 ' 

2 d\ny 

We define the temperature bins (T, T + AT) so that there are at least two clusters in each. 
For a fixed background cosmology and power law index a, the left-hand side of equation ( p3|) 
contains only observed quantities, and the right-hand side is a function of y and the PDF. 
For each PDF, we therefore can solve for y implicitly as a function of T. We then convert T 
to R using equation (|23|) and our derived value of y, and obtain a value of cr^ ustcrs = S c /y as 
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a function of the value of R inferred from each temperature bin. Multiplying by ^q 6 gives 
the cluster-inferred ^^ usters ; which can be compared directly with the values from the rms 
measurements. 

A clear comparison with rms-inferred 77 can be made if we extrapolate all values of r\ 
to a common scale, using the power law model we have assumed throughout. We choose 
10ft. _1 Mpc because this scale corresponds to a temperature within the range of the obser- 
vations for Qq = 0.3 — 1.0: 8.6 keV for Q Q = 1 and 3.4 ~ 4.0 keV (depending on <j)) for 
Qo = 0.3. For Qq — 0.3, the scale R = 8/i _1 Mpc corresponds to a temperature of 2.2 ~ 2.6 
keV, outside the range of observations we consider. 

For each temperature bin, we therefore extrapolate the derived value of rjn to 10/i _1 Mpc 
using r]f> oc or qc R~ a , and obtain a measurement of 7710 for each bin. We wish to find the 
resulting probability distribution in r/io- We find that the statistical errors in 7/10 from equa- 
tions (^) and (|33"D are smaller than the scatter of the values of r/10, though the discrepancies 
are not as drastic as in the case of the rms measurements. Certainly, then, there remains 
some type of systematic error. We sum the distribution functions of each extrapolated r)i 
to obtain the total distribution function, assuming each error is Gaussian: 

P(77 10 |£ c i us PDF) = -£^=e~H - J , (34) 

where the sum is over the temperature bins labeled by i. Note that we are summing over 
the individual distributions because the individual temperature bins are not independent 
measures of r/10; i.e., the data set as a whole is taken as a single measurement, with the error 
given by the total distribution "P(r/io|-Dcius PDF). 

Figures |5|-|7| show the cluster-inferred 7710 distributions for power laws a = 0.83, 0.60, and 
1.10, and for different values of Qq and the PDF. Also shown are the likelihoods for mean 
peculiar velocity/redshift distortion-derived rms values discussed in §2, and the Willick et 
al. (1997b) result alone. The means and standard deviations of the resulting distribution of 
7710 measurements are given in Table 0, for various values of Qq and a. 

The most prominent feature of the inferred rms from clusters is its strong dependence on 
the PDF. Gaussian and texture models give clearly different predictions for any given cluster 
data set. The degree to which the cluster-inferred value of 7/10 differs between Gaussian and 
texture PDFs is directly related to the degree to which the tail of the texture PDF is 
distinguishable from a Gaussian. 



4. Discussion 
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4.1. Which PDF is More Likely? 



We quantify the comparison between the rms and cluster-inferred determinations of t]iq 
with a Bayesian model comparison analysis (Loredo 1990). For now, we fix the background 
cosmology and the power law index a, and suppress the dependence on these assumptions. 
We ultimately seek the probability that a certain PDF is correct, given the cluster data -D c ius 
and the rms data D rms : 

P(PDF) 



V(PDF\D clus A 



-V(D clus LUJPDF), 



(35) 



^(^clus ^rms 1 

where we have used Bayes Theorem. The fraction in front consists of a prior (which we 
take as uniform) and a normalization factor. Using the product rule, we separate the last 
probability: 

V(D clus -D r ms I PDF) = V(D clus \PDF)V(D rms \D clus PDF). (36) 

Since we are primarily interested in the rms-versus-tail comparison, we take the first proba- 
bility in equation (|36|) to be independent of the PDF. A more sophisticated analysis would 
evaluate the "goodness of fit" of the cluster data for each PDF, and thereby estimate the 
absolute probability P(D clus |PDF). However, as mentioned in §3, the scatter in the inferred 
values of ?7io, c ius is greater than the calculated Poisson errors, so any such analysis will have 
to take into account unknown systematic errors (e.g. using the model described in §2.2). 



The dependence of D ims on -D cms in the second probability in equation ( |3T)| ) is through 
the value of 7] w ; using the summation rule gives 



V(D imB \D clus PDF) = J dr) W V(D ims \r) 10 PDF)V(r) 10 \D clus PDF). 



(37) 

The resulting expression for equation (|35|), suppressing all factors assumed to be independent 
of the PDF, is 



(38) 



V(PDF\D clus D rms ) cc J d V ioV(D ims \T] 10 PDF)V(r]w\D clns PDF). 

The cluster distribution function in r] W , V(r]i Q \D c \ ns PDF), is given by equation (|34T). 
The rms probability depends only the power law slope a, and the weighting of the different 
analyses, and is independent of the PDF. If we use the error model in §2.2, the probability 
density P(-D rms |?7 10 PDF) is given by the likelihood function from equation ([□]) and Figure B 

V(D ims \ Vl0 PDF) = C ims (Vio)- (39) 

If we are considering only the Willick et al. (1997b) data, -D-wniick; and we assume normal 
errors, we use 



7>(Ayiiiick|r/io PDF) 



1 



2ne 2 w 



exp 



1 f Vio -Vw \ 2 

2 V e w / 



(40) 
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where r^w and ew are given by the values inferred from their analysis, also listed in Table |I[ 

Since we do not know the full range of possible PDFs, we take the Bayesian ap- 
proach of considering the likelihood ratio (the "odds", or "Bayes factor" for uniform pri- 
ors, Loredo 1990) of a Gaussian versus a texture PDF. The values of the likelihood ratio 
■p(Gaussian)/P(texture) are listed in Table H, for various background cosmologies and val- 
ues of a. If we take the average of the rms measurements, a Gaussian PDF is favored for 
any background cosmology. The probability of a Gaussian PDF is always > 2x that of 
a texture PDF. The confidence level of the accepting a Gaussian hypothesis is 70 — 88%, 
roughly la ~ 1.5a. If we only use the Willick et al. (1997b) measurement, however, then 
which PDF is favored depends strongly on the background cosmology, and somewhat on the 
power law index a. In the case of a = 0.83, if we live in an open Qq = 0.3 universe, then the 
Gaussian PDF is 3x more likely than is the texture PDF; if we have a flat Qq = 0.3 universe 
(A = 0.7), then a Gaussian and texture PDF are roughly equally likely; and if Qq = 1, a 
texture PDF is 3x more likely than a Gaussian. 



4.2. Constraints from Cluster Evolution? 

We have only considered the present day abundances of clusters, and we have a degen- 
eracy between the assumed PDF and the background cosmology. The evolution of cluster 
abundances, however, offers an opportunity to break this degeneracy. In particular, high Qq 
universes evolve more rapidly than those with low Qq. Gaussian models with Qq — 1 are 
already heavily disfavored because they predict too rapid an evolution in cluster abundances 
out to z> 0.3 (Henry 1997; Donahue et al. 1997). The fact that the Willick et al. (1997b) 
measurement leads to the Qq = 1 Gaussian scenario being disfavored relative to a texture 
scenario simply reinforces the case against a flat matter-dominated universe for Gaussian 
fluctuations. 

The formalism of §3 relating cluster abundances to the PDF and the value of r\ is 
fully generalizable to z > 0. As an example, we show in Figure [8] the cumulative number 
density of clusters > 10 keV, as a function of redshift for various scenarios. These curves are 
normalized to a present-day number density of 1 x 10 -8 /i 3 Mpc~ 3 (see Figure ^). For fixed 
background cosmology, the texture scenarios evolve slower than the Gaussian ones because 
of the shallower texture PDF. For fixed PDF, the evolution of the number density evolves 
fastest with Qq = 1 and slowest with open Qq = 0.3. A Gaussian PDF with a flat Qq = 0.3 
background cosmology has roughly the same evolutionary history as a texture PDF with 
Qo = 1. If evolution is observed to be significantly slower or more rapid than this, then an 
Q = 1 texture scenario can be ruled out. 
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Henry (1997) reports that the number density of z — 0.3 clusters at the hot end of 
the temperature function ( £ 6.5 keV) is a factor of ~ 2.5 lower than today (his Figure 2), 
though this is still within the calculated scatter due to Poisson statistics. For Qq = 1, the 
Gaussian model predicts a factor ~ 8 decrease in the density of such clusters at z — 0.3 
while the texture model predicts a factor ~ 3 decrease. For the Qq = 0.3 scenarios, the 
reduction factor is 1.1 ~ 2. Henry (1997) reports that the cluster evolutionary data rule out 
the Gaussian f2 = 1 at 99% confidence. 

Donahue et al. (1997) report an extremely hot X-ray temperature of 14.7 ± 4 keV for 
cluster MS1054-03, at redshift 0.828. They indicate that this result, along with data from 
other high-redshift clusters they have studied, implies that, within errors of ~ 50%, the 
number density of clusters with T > 10 keV is essentially unchanged from z = to z = 0.8. 
If this result is confirmed, even a texture PDF would be ruled out for fi = 1, since it 
predicts a decrease in the number density by a factor of ~ 30. Further observations at 
redshift z = 0.5 ~ 1.0 will almost certainly resolve the ambiguity between models with 
different Qq. 



5. Summary and Conclusions 

We show that comparing the amplitude of tail fluctuations indicated by the abundances 
of clusters of galaxies, to the amplitude of the rms mass fluctuation as indicated by galaxy 
peculiar velocities and redshift space distortions, can be a strong test of the Gaussianity of 
the PDF of the initial density field, with a weak dependence on the slope of the density 
fluctuation spectrum. 

Observations of peculiar velocities and redshift distortions in the linear regime can 
measure the rms level of gravitational potential fluctuations present in the universe in the 
combination i]r = <7rQq S , independent of galaxy bias. Various existing measurements of 
this quantity are not consistent with one another within their stated errors. We assume 
the existence of a systematic error in each measurement, and use a maximum likelihood 
technique to combine these recent measurements of The systematic error leads to a 
~ 20% uncertainty in the value of i] R . This procedure gives r]i rms = 0.39±0.07 — 0.51±0.10, 
depending on the assumed slope of the mass fluctuation spectrum. We use this value, and 
one recent measurement by Willick et al. (1997b), which indicates r] W = 0.27 — 0.30 ± 
0.04, to compare with the cluster-inferred measures of 7710- We chose the Willick et al. 
(1997b) measurement because it is the most recent analysis which explicitly filters close to 
the 10/i~ 1 Mpc scale on which clusters form. 
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The value of r] 10 inferred from cluster abundances is calculated using the Press-Schechter 
approach for both a Gaussian and a texture PDF. We find a simple relation for the typical 
redshift of formation of clusters in both models. This is particularly important in the texture 
scenario, where the shallower tail of the PDF leads to significant early formation of clusters. 
This typical redshift is important in the relation between the virial temperature of the 
clusters and the comoving radius of the initial perturbation. For a texture PDF, the observed 
abundances of X-ray clusters implies r] w = 0.20 ± 0.03 — 0.25 ± 0.04 for Q = 0.3, and 
7]iq = 0.22 — 0.29 ± 0.03 for Qq = 1.0. Given a Gaussian PDF, on the other hand, the cluster 
abundances indicate r] 1Q = 0.29 ± 0.04 - 0.33 ± 0.05 for Q = 0.3 and r] w = 0.34 ± 0.05 - 
0.42 ± 0.04 for Q = 1. 

We calculate the relative likelihoods of each PDF for a range of cosmological parameters. 
Using the maximum likelihood average of the rms measurements with the cluster abundances 
implies the Gaussian model is favored for any Qq. However, using the Willick et al. (1997b) 
rms measurement with the clusters implies that a Gaussian is favored for open Qq = 0.3, a 
Gaussian and a texture PDF of roughly equal likelihood for flat Q Q = 0.3, and a texture PDF 
is favored for fl = 1. In each case, the accepted hypothesis is accepted at a ~ 1.5<r level, 
with likelihood ratio £ 3. Preliminary constraints from cluster evolution indicate that low 
Qo models are probably more viable than f2o=l models. Our conclusions can are summarized 
in Table |], where a \J (X) indicates the PDF is favored (disfavored) over the alternative by 
3x. A ? next to the ^/or X indicates the factor is ^3. 

Our results imply that the textures+CDM scenario faces some challenges from obser- 
vations on 5 — 25/i -1 Mpc scales, completely independent of CMB measurements. Although 
the current state of the measurements of the rms prevents us from firmly distinguishing the 
Gaussian and texture models, the principle behind this method is clear. The observed cluster 
abundance at z = predicts significantly different values for 7] w for the two models. More 
accurate direct determinations of the measured rms from future redshift surveys and peculiar 
velocity measurements, or from an improved understanding of the systematic effects which 
are present in the current analyses, will surely enable us to definitively test any distribution 
function of initial density perturbations. Furthermore, upcoming high redshift X-ray cluster 
data will likely break any remaining degeneracy due to the assumed background cosmology. 
Soon, perhaps, we will finally be able to distinguish different models for the initial PDF of 
the universe at 10/i _1 Mpc. 

We thank Avishai Dekel and Ed Turner for useful conversations. MAS acknowledges the 
support of the Alfred P. Sloan Foundation, Research Corporation, and NSF Grant AST96- 
16901. 
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A. What is the Correct Correspondence Between Tophat and Gaussian 

Filtering? 

One of the traditional ways of relating different filtering schemes is by matching the 
volume over which they filter. For tophat and Gaussian, then, this would correspond to 
AttR^/3 = (2ttRq) 3 ^ 2 , or R T R g /1.Q. While this might be appropriate when relating 
the mass of objects picked out by these filter scales, it is not the best way of relating the 
wavelengths of fluctuations measured by each filter. Here we derive a method of relating 
tophat- and Gaussian-filtered fields by their maximum correlation. 

Consider a homogeneous, isotropic scalar field -F(x) and its Fourier transform in the 
following convention: 

F(x) = ^F(k)e^ = j^jd 3 kF(k)^- (Al) 

F(k) = 1 J d 3 xF{x)e~ ik - x . (A2) 

Let the power spectrum P(k) = (|F(k)| 2 ). The (dimensionless) tophat and Gaussian filters 
and their Fourier transforms are then 

3 

W Rt (x) = V — -3 if |x| < R T , otherwise, (A3) 

W RT (k) = 33 f D RT) = W T (kR T ), (A4) 
kRt 

Wrg(x) = V (2nR* G )^ eM ~ x2/{2R2G)) ' (A5) 

W RG (k) = exp(-k 2 R 2 G /2) = W G (kR G ), (A6) 

where ji(x) is the spherical Bessel function of order I. The filtered fields are F Rtg (x) = 
V^ 1 J d 3 yF(x. — y)WR T G (y). Using the convolution theorem, the variance for tophat and 
Gaussian filters is: 

*t,g( R t,g) = J^J dk ^k 2 P(k)Wl G (kR T , G ). (A7) 

How do we relate results with a tophat filter and those with a Gaussian filter? We maxi- 
mize the correlation between the two filtered fields, normalized by their respective standard 
deviations: 

(A8) 
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For Gaussian Random Fields, the one-point distribution of each of the filtered fields is 
Gaussian, and the joint distribution of the two fields is a bivariate Gaussian. In this case the 
quantity p{Rr, Rg) is precisely the normalized correlation coefficient of the joint distribution 
of the two fields. The two filtering schemes, therefore, are most correlated when p(Rt, Rg) is 



maximized. In Table [Al], we list the peak correlations for power-law power spectra P(k) oc k n 
for n = —2,-1,0, and 1. For decreasing power spectra (n < 0), the correlation between 
tophat and Gaussian filtering is very good, greater than 87%, with the ratio Rt/Rg — 1-75 ~ 
2. 

It should be noted that we are not comparing the value of a between the two filtering 
schemes — we are seeking the scales at which the two filters are most correlated. For example, 
for a CDM, Qq = 1, n — 1, power spectrum, using the maximum correlation, 12/i _1 Mpc 
Gaussian smoothing corresponds to 22.7/z _1 Mpc tophat smoothing; but ctg(12)/ctt(22.7) = 
0.86, while a G (12)/<r T (25.4) = 1.0. 
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Fig. 1.— RMS mass fluctuations as inferred from peculiar velocity and redshift dis- 
tortion analyses, as described in §2, and listed in Table ffl. The solid line is T]r = 
0.44(i2/10.0/i~ 1 Mpc) -0 * 83 , where a is that for IRAS galaxies, and 7710 is taken from the 
maximum likelihood average in Table |I[ 
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Fig. 2. — Probability density for the error analysis in §2 of the rms measurements of 7710, for 
power law indices a =0.60, 0.83, and 1.10. The thin dotted line is a Gaussian distribution 
with the same mean and 68% confidence interval (la error) as the probability density for 
a = 0.83. 
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Fig. 3. — Tails of the effective distribution functions fP(y), for a Gaussian (solid) and 
texture (dashed) PDF. 
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Fig. 4. — Cumulative temperature function from Henry & Arnaud (1991) (thick solid line), 
and from the XBACS survey QEbeling et al. 19961) (thick dashed line) , both with la Poisson 
errors from equation (|32|) (thin histograms). 
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Fig. 5. — Cluster-inferred 7710 compared with peculiar velocity/redshift distortion-inferred 
7710 at 10/i _1 Mpc for a = 0.83, for various assumed background cosmologies. The solid-filled 
distribution functions are for a Gaussian PDF, and the dashed-filled distribution functions 
are for a texture PDF. The filled triangle is the maximum likelihood average of the measure- 
ments in Table [l]. The open pentagon is the Willick et al. (1997b) result. 
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Fig. 6. — Cluster-inferred i] w compared with peculiar velocity/redshift distortion-inferred 
7710 at 10/i _1 Mpc for a = 0.60. The solid-filled distribution functions are for a Gaussian 
PDF, and the dashed-filled distribution functions are for a texture PDF. 
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Fig. 7. — Cluster-inferred i] w compared with peculiar velocity/redshift distortion- inferred 
7710 at 10/i _1 Mpc for a = 1.10. The solid-filled distribution functions are for a Gaussian 
PDF, and the dashed-filled distribution functions are for a texture PDF. 




Fig. 8. — Evolution of the number density of clusters with temperature > 10 keV, normal- 
ized to the present day abundance, for Gaussian and texture PDFs, and various background 
cosmologies. 
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Table 1. RMS Measurements of rjio 



Reference 


R 


Vr 




mo 










a = 0.60 


a = 0.83 


a = 1.10 



Zaroubi et al. 1997 



IKoJatt & Dekel 19961 



WilJick et al. 1997b! 



Sigad et al. 1997 



Fisher et al. 1994b 



Cole et al. 1995a,b 
Fisher, Scharf, fc Lahav 1994 
Maximum Likelihood Av 



10 

23 

8 

23 
12.5 
35 
30 
10 



0.70 ±0.09 
0.32 ±0.11 
0.34 ±0.05 
0.25 ±0.04 
0.21 ±0.13 
0.11 ±0.03 
0.20 ±0.03 



0.70 ±0.09 
0.53 ±0.18 
0.30 ±0.04 
0.41 ±0.07 
0.25 ±0.15 
0.22 ±0.06 
0.39 ±0.06 
0.39 ±0.12 



0.70 ±0.09 
0.64 ±0.22 
0.28 ±0.04 
0.50 ±0.08 
0.26 ±0.16 
0.30 ±0.08 
0.50 ±0.07 
0.44 ±0.12 



0.70 ±0.09 
0.80 ±0.27 
0.27 ±0.04 
0.62 ±0.10 
0.27 ±0.17 
0.42 ±0.12 
0.67 ±0.10 
0.51 ±0.15 



Table 2. XBACS Cluster Measurements of i]i 



Cluster S 


cenario 




^10 








a = 0.60 


a = 0.83 


a = 1.10 


Gaussian fl 


= 1 


0.423 ± 0.036 


0.381 ±0.041 


0.342 ±0.049 


texture Qo = 


1 


0.287 ±0.033 


0.252 ±0.031 


0.220 ±0.032 


Gaussian fl 


= 0.3 fiat 


0.326 ±0.035 


0.325 ±0.038 


0.329 ±0.047 


texture Q = 


0.3 flat 


0.245 ±0.037 


0.235 ±0.034 


0.231 ±0.036 


Gaussian fl 


= 0.3 open 


0.298 ±0.029 


0.293 ±0.033 


0.293 ±0.042 


texture Qq = 


0.3 open 


0.215 ±0.028 


0.204 ± 0.027 


0.198 ±0.029 
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Table 3. Probability of Gaussian (G) Relative to Texture (T) PDF 



rms 


fl and 




P(G)/P(T) for a = 




Averaging 


Geometry 


0.60 


0.83 


1.10 


Maximum 


1.0 


2.28 


6.52 


7.00 


Likelihood 


0.3 flat 


3.28 


4.18 


4.41 




0.3 open 


3.98 


4.71 


4.62 


Willick et al. 


1.0 


0.124 


0.318 


0.704 


(1997b) only 


0.3 flat 


1.53 


1.25 


0.789 




0.3 open 


3.28 


3.01 


2.13 



Table 4. Summary of Tests of Gaussian (G) Versus Textures (T) 



Tests 


^0 


= 1 


^0 = 


0.3 flat 


^0 


= 0.3 open 




G 


T 


G 


T 


G 


T 


rms average and z = clusters 


V 


X 


V 


X 


V 


X 


Willick et al. (1997b) and z = clusters 


X 


V 


?— 


?— 




X 


cluster evolution 


X 


?x 


V 


V 


V 


V 



Note. — A u y/" ("X") indicates the PDF is favored (disfavored) over the alternative by 
^,3x or £ 1.5a; a "?=" indicates roughly equal likelihood; a "?" next to the "-^/" ("X") 
indicates a ^ 3 or <, 1.5<r result in favor (against). 
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Table Al. Correspondence Between Tophat and Gaussian Filtering 



index n 


Rt/Rg 


p(R T , Rg) 


(P(k) oc k n ) 


(maximum correlation) 




-2 


1.94 


0.99 


-1 


1.85 


0.97 





1.75 


0.87 


1 


1.66 


0.57 



